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1. 


Introduction 


Spherical  Diffusion  Flames 


Earlier  chapters  have  been  concerned  with  flames  for  which  the  reactants  are 
supplied  already  mixed.  When  the  reactants  are  intially  separated  and  diffuse  into 
each  other  to  form  a  combustible  mixture,  the  ffrljme  is  called  a  diffusion  flame. 

A  bunsen  burner  with  its  air  hole  closed  supports  a  diffusion  flame  between  the 
gas  supplied  through  the  tube  and  the  surrounding  oxygen-rich  atmosphere.  A  ? 
candle  supports  a  vapor  diffusion  flame,  so-called  because  the  fuel  is  produced, 
through  liquefaction  and  subsequent  evaporation  of  the  wax,  by  the  heat  of  the  flame 
itself. 

I 

Premixed  and  diffusion  flames  share  certain  features,  but  the  differences 
are  more  striking  than  the  similarities.  For  example,  there  is  no  unlimited  plane 
flame,  as  can  be  seen  from  Sec.  II. 7.  The  oxidant  fraction  would  be  constant 
beyond  the  flame  sheet,  at  which  the  value  would  be  smaller  than  upstream.  But 
a  diffusion  flame  requires  the  oxidant  fraction  to  be  zero  upstream,  since  other¬ 
wise  there  is  premixing,  so  that  a  contradiction  is  reached.  A  similar  argument 
applies  when  the  fuel  is  supplied  at  a  finite  point,  where  now  it  is  the  oxidant 
flux  which  must  be  zero.  Moreover,  since  cylindrical  flames  are  just  geometrically 
attenuated  versions  of  plane  flames  (Ch.  VII),  there  is  no  cylindrical  diffusion 
flame  either  (Ludford  1977). 

In  order  to  gain  some  insight  into  the  nature  of  diffusion  flames  we  seek 
a  simple  one-dimensional  configuration  which  bears  some  relation  to  physical 
reality.  One  possibility  is  to  introduce  the  oxidant  at  a  finite  point  of  the 
plane  flame  but,  curiously  enough,  details  of  that  have  never  been  worked  out. 

To  be  sure,  the  chambered  diffusion  flame,  in  which  the  fuel  supply  is  also  at 
a  finite  point,  has  recently  been-  considered  by  Matalon,  Ludford  and  Buckmaster 


(1978).  But  that  is  already  too  complicated  for  our  purposes. 

Two  choices  have  been  popular  in  the  literature.  One  is  the  counterflow 
diffusion  flame  (Fig.  l)  in  which  opposing  jets  of  fuel  and  oxidant  collide, 
supporting  a  locally  one— dimensional  flame  in  the  neighborhood  of  the  stagnation 
point.  The  configuration  has  the  advantage  of  being  relatively  easy  to  set  up 
experimentally  but  suffers  analytically  from  the  gross  compromises  which  have  to 
be  made  with  the  fluid  mechanics.  For  that  reason,  and  because  of  its  similarity 
with  the  alternative  to  be  discussed  next,  we  shall  not  consider  the  counterflow 
diffusion  flame. 

The  second  choice  is  the  spherical  diffusion  flame  which,  because  it  occurs 
in  the  burning  of  fuel  drops  (Fig.  2),  has  great  technological  interest.  As  for 
the  candle,  it  is  then  a  vapor  flame,  the  liquid  fuel  being  vaporized  at  its 
assumed  spherical  surface  by  heat  conducted  back  from  the  surrounding  flame  where 
the  vapor  burns  with  the  ambient  oxidant.  Although  the  size  of  the  drop  continually 
decreases  as  fuel  is  consumed,  the  quasi -steady  approximation  (Williams,  1965) 
neglects  this  effect.  Such  steady  burning  can  be  produced  experimentally  by 
forcing  the  fuel  through  a  porous  sphere  in  a  gravity- free  environment.  Indeed, 
by  supplying  the  fuel  already  as  a  vapor  the  whole  range  of  spherical  flames 
can  be  examined. 

If  these  phenomena  are  studied  for  the  express  purpose  of  solving  the 
practical  problem  from  which  they  derive,  then  one  has  to  be  concerned  with  the 
validity  of  the  modelling.  Thus  many  studies  of  the  fuel-drop  problem  have 
'examined  amongst  other  things  the  validity  of  the  quasi-steady  approximation, 
about  which  reservations  have  been  expressed.  Such  practical  considerations 
are  not  our  concern,  since  we  regard  the  models  simply  as  mathematical  idealizations, 
albeit  with  physical  reality  continually  in  mind,  whose  study  can  provide  some 
insight  into  the  nature  of  diffusion  flames.  The  experimentalist  may  then  feel 


challenged  to  devise  real  flames  which  model  the  mathematics,  as  in  the  porous-sphere 
experiment  mentioned  above. 


We  have  to  deal  with  the  spherically  symmetric  form  of  the  equations  ( I .  i+l-l+I; ) , 
in  particular 
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where  the  radial  velocity  v  is  the  only  component,  T  =  1/p,  t=  t/0 
and  (if  the  radius,  a,  of  the  supply  sphere  is  used  to  define  units) 

(4)  n  =  D  Y^Y^2  e~9/T 

with  D  given  by  the  formula  (1.1+7).  The  equations  for  the  product  and  inert  species, 
which  have  not  been  written,  are  complicated  by  the  implicit  assumption  that 
is  not  necessarily  equal  to  *£g.  ■^ie  mom®iitum  equation,  which  has  also  been  omitted, 
plays  no  role  other  than  determining  the  pressure  field  when  there  is  spherical 
symmetry.  Throughout  this  chapter  we  shall  take 

(5)  V1  *  v2  =  1  80  that  °l  =  °2  = 

(corresponding  to  a  bimolecular  reaction)  for  the  sake  of  simplicity,  even  though 
retaining  general  values  would  cover  the  corresponding  analysis  in  Ch.  VII 
*  1,  v2  *  0)  showing,  in  particular,  why  there  are  so  many  similarities. 

Finally,  unsteady  terms  have  been  retained  with  a  view  to  stability  considerations 
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on  the  long  time  t.  To  avoid  subscripts  we  shall  change  the  notation  to 

(6)  Y  *  Y,  Yg  =  Z,  Jx  =  J,  J2  =  Ki^±  ~Xt  X2  =~ffl  . 

All  six  values  T_.  Y_,  Z„,  Tl,  Y',  Z*  may  be  assigned  at  the  supply  sphere, 
s  s  s  s  s  s 

In  varying  the  mass  flux  M  from  the  sphere  we  shall  in  fact  suppose  that 
M~^T' ,  M-^Y' .  M-^Z*  are  prescribed,  these  being  fluxes  per  unit  mass  of  mixture 

8  3  S 

supplied.  (The  same  was  true  for  the  plane  flame;  M  was  then  hidden  in  the 
coordinate.)  It  is  more  convenient  to  suppose  that  the  two  fluxes 

(7)  J,  -  I,  -  .  K3  =  Za  -  M -\Mtl  ' 

are  given.  Clearly  a  pure  diffusion  flame  requires  J  >0,  K  =  0  so  that  the 

s  s 

sphere  is  a  source  of  fuel  but  not  of  oxidant;  for  simplicity  again,  we  shall 
take. 

t8)  Ja  =  1,  Ks  =  0 

making  the  sphere  purely  a  source  of  fuel.  Likewise  the  values  T  ,Z  at  infinity 

00  00 

will  be  given,  leaving  only  T  ,M  ^T'  =  L  among  the  six  original  parameters. 

s  s 

In  addition  Y  may  be  assigned,  and  we  shall  take 
00 

(9)  I.  -  0 

to  ensure  that  all  the  fuel  has  originated  at  the  supply.  There  are  then  seven 
boundary  conditions  on  the  sixth-order  system,  so  that  we  may  expect  M  to  be 
determined  by  D. 

2.  Steady  Combustion  for  ^  =  1.  DamkOhler  Humber  Asymptotics. 

The  qualitative  nature  of  the  steady  solution  is  not  affected  by  the  Lewis 
numbers;  it  is  only  in  the  stability  question  that  they  play  an  important  role. 


A 


To  simplify  the  discussion,  which  will  be  complicated  enough,  we  shall  therefore 
assume  for  the  present  that  both  Lewis  numbers  are  unity. 

In  the  steady  state  the  continuity  equation  can  be  integrated  to  give 


(10) 


o  2 

r  pv  =  M  (const.)  or  pv  =  M/r  . 


krr&^M  is  the  rate  at  which  fuel  is  supplied,  but  is  is  only  the  burning  rate 

if  the  total  outflow  of  fuel  at  infinity  is  zero.  Nevertheless  M  is  loosely  called 

the  burning  rate,  and  we  shall  adopt  the  term  also. 

The  Shvab-Zeldovich  variables  Y^  =  T-Y^/a^  satisfy  f 


(ID 


M  _d  ] 
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=  0 


for  1<  r<  00  . 


r  r 

Since  all  solutions  are  bounded  at  infinity,  the  analysis  immediately  differs  from 
that  for  the  plane  (premixed)  flame.  As  already  anticipated,  all  six  basic 
parameters  TQ,  L,  JQ,  KQ,  T^,  ^»may  be  assigned  to  give,  in  the  case  (8) 


(12)  T+2Y  =  +  (Ta-Tj(l-e“M/r), 

(.13)  T+2Z  =  (Ts-L)(l-e“M/r)  +  (T  +  2Y^)e-M/r 


where 


(14)  Ta  =  Ts-L+2 


is  the  adiabatic  flame  temperature .  The  problem  now  reduces  to  one  for  T  alone , 

namely 

^(T)  *  DYZ  e_0/T  for  1  <  r  <  » 

T  ,  IijT^  prescribed, 
s 


(15) 
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Since  there  are  three  boundary  conditions  for  this  second-order  differential 
equation,  M  is  determined  as  a  function  of  D  as  part  of  the  solution. 

Until  the  advent  of  activation-energy  asymptotics,  analytical  discussion 
was  preoccupied  with  Damkohler-number  asymptotics  (i.e.  D->  0  or  <*») .  Knowledge 

of  these  limits  is  of  some  importance,  so  they  will  be  briefly  discussed  first. 

Frozen  combustion  occurs  when  the  chemical  reaction  is  very  weak,  i.e. 

D  -*-0.  The  heat  conducted  back  to  the  supply  sphere  then  comes  from  the  resevoir 
at  infinity  (provided  T#  >Ts)  and  the  limit  solution  is 

(16)  T  =  (Ts-L)  +  LeM(l'l/r)  f 

where 

0-7)  M  =  =  m  [1  +  (T^  -  Ts)/L  ] 

We  shall  only  be  concerned  with  conductive  heat  sinks,  i.e.  L>  0,  so  that 
CIS)  Tg<  T, 

is  required  for  Y  to  be  positive;  then  M  is  also  positive,  as  expected. 

The  result  (16)  is  not  uniformly  valid  since  for  large  values  of  r  the 
neglected  reaction  term  decays  more  slowly  than  the  transport  terms.  Defining 

2  -0/T» 

(19)  e  =  Z  e 

00 

leads  to  the  expansion 

(20)  T  =  T.-  [T^-  Ts  +  L  -  2(l-e"R)]+  o(e) 

for  the  coordinate  R  =  er.  The  temperature  gradient  far  from  the  supply  is 
thereby  corrected,  so  that  the  heat  flow  from  the  infinite  reservoir  can  be 


calculated  as 


(21)  lim  (4m-2  f£)  =  U,M(T  -T  ). 

__  dr  00  a. 

When;  "** 

>T^  an  0(1)  amount  of  heat  is  supplied  to  the  environment,  which  is  usually 
the  goal  of  combustion.  The  frozen  limit  is  not  an  extinguished  state,  but  rather 
one  in  which  all  the  reaction  takes  place  at  essentially  constant  temperature  far 
from  the  supply. 

The  equilibrium  limits  are  obtained  as  D  ■*“.  They  are  of  great  practical 
interest  since  even  under  standard  atmospheric  conditions  many  diffusion  flames 
do  have  large  Damkohler  numbers.  Clearly  the  limit  D-*”  is  singular,  since  tne 
highest  derivative  in  equation  (15)  is  multiplied  by  a  vanishingly  small  parameter 
D-^.  We  may  anticipate  thin  regions  where  the  second  derivative  is  very  large 
and,  when  located  in  the  interior  of  the  combustion  field,  they  are  called  Burke- 
Schumann  flame  sheets  (Kassoy  &  Williams  1968).  They  should  not  be  confused 
with  the  flame  sheets  that  occur  in  activation-energy  asymptotics.  The  thin  regions 
can  also  be  boundary  layers  but  we  shall  look  at  the  other,  more  common  possibility 
first. 

Outside  the  flame  sheet  the  limit  equation  is  simply 


(22) 


YZ  =  0, 


which  represents  chemical  equilibrium  for  the  irreversible  reaction.  The 
combustion  field  is  divided  into  regions  where  either  Y  or  Z  vanishes.  Since 
fuel  is  supplied  at  the  sphere  and  oxidant  at  infinity  it  is  natural  to  put  an 
oxidant-free  region  contiguous  with  the  sphere,  separated  at  r  =  r#  by  a  flame 
sheet  from  a  fuel-free  region  extending  to  infinity.  The  Shvab-Zeldovich  relations 
(12,13)  then  imply 
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f  *1 

(23) 

Z  = 

j 

(24) 

Y  = 

and  the  temperature  at  the  sphere  is  T_  only  if 


r  T  -  T  +2Z 

(25)  M  =  Me  =  SLn\l  +  -^-S - —  J  . 

Continuity  of  temperature  at  the  flame  sheet  shows  that  it  must  he  located  at 

(26)  r*  =  M Jin  (l  +  Z  ), 

“  00 

so  that  the  flame  temperature  is 

(27)  T*  =  T«  +  ZjTa-Tj(l  +  Zj. 

These  equations  are  the  essence  of  the  so-called  Burke-Schumann  solution, 
although  it  is  still  necessary  to  show  that  there  is  a  structure  linking  the  two 
sides  of  the  flame  sheet.  We  shall  he  content  with  setting  up  an  apparently 
well-posed  problem  for  this  structure. 

It  is  clear  from  the  oxidant  species  equation  (3b)  that,  for  D  large  hut 
not  infinite,  there  can  be  no  algebraic  perturbation  of  Z  =  0  in  r  <  r#.  The 
temperature  is  therefore  given  to  all  orders  by  the  formula  (23)  so  that,  on 
writing 

_ T  /o  *1  /-n  —  0/T|£ 

(28)  M  =  MQ  +  D,  Mx  +  o(d;XM),  Mq  =  Me,  D*  =  ?gDe 

(a  choice  justified  by  the  consistency  of  the  structure  problem  below),  we  find 

(29)  T  =  T~  +  D;1/3M1(Ts-Too-L-2Zoo)e  ^  //V>(d;1/3)  . 

where  T”  is  the  original  T  in  r  <  r#  with  M  replaced  by  MQ.  Similarly, 
there  is  no  algebraic  perturbation  of  Y  in  r  >  r#  so  that  the  temperature  is 


r*  =  M Jia  (1  +  Z  ), 
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-M  Jr 


(31) 


T  =  T*  +  D“1//3M1(Ta-T  Je  °  /r  +  o(D#^3) . 


Within  the  flame  sheet  we  set 
(32)  T  =  T#  +  D;l/3t  +  o(d;1/3),  r  =  r,  +  D;l/3£ 

to  obtain 


(33) 


“W  =  -  -  £)]  It  ♦  (j^=-  a)(-|-  -  |)1  . 

•  r#  00  £ 


Matching  with  the  expansions  outside  requires  that 

T_-T_  M05  M1 


(34) 

(35) 


t  =  (- 


oo  a 

I+Z 


+  2) (— r—  -  — )  +  o(l)  as  £-»■  -«  , 


T  -T  Mn 

00  a  u  .  /  •%  \  _  . 

U  "  1+Z  ZS  0<1)  “ 


dt 


oo  r; 


It  is  reasonable  to  expect  that,  for  some  unique  in  the  initial  condition 

(34)  ,  integration  of  the  structure  equation  (33)  from  g  =  will  lead  to 

the  correct  slope  (35)  at  £=  +<=°  .  That  indeed  is  found  numerically,  which  is 
the  only  way  of  carrying  out  the  integration.  We  may  therefore  assert  that 
there  is  a  structure  and  that  its  computation  gives  the  perturbation  of  the 
burning  rate. 

The  Burke-Schumann  solution  does  not  hold  for  all  parameter  values.  An 
underlying  assumption  is  that  a  flame  sheet  forms  in  the  interior  of  the 
combustion  field,  so  that  for  consistency  r#  must  be  greater  than  1,  i.e. 


(36)  V  Ts>  (L-2)Zw  . 

Then  MQ  is  automatically  positive,  as  might  be  expected.  As  equality  is 
approached  the  flame  sheet  moves  to  the  surface  r  =  1  and  our  analysis  breaks 
down.  Although  the  resulting  surface  flame  (Buckmaster  1975)  does  not  play  as 
important  a  role  in  practice  as  the  Burke-Schumann  flame,  a  brief  description  of 
it  will  be  given  here. 
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When  the  flame  sheet  is  a  boundary  layer  the  oxidant-free  region  is  absent 
and  only  the  fuel-free  region  (24)  remains.  Since  the  temperature  must  tend  to 
T  as  r  -*-1, 

3 

(37)  M  =  Mg  =  toll  +  (Tw-  T?)/vL-2)]  , 


a  result  quite  different  from  the  Burke-Schumann  value  (25)  but  to  which  it  is  eqaui 
at  the  inequality  (36).  The  result  (37)  is  that  for  weak  burning  with  L  replaced 
by  L-2,  where  the  2  can  be  identified  as  the  heat  released  in  the  flame. 

In  addition,  the  flame-sheet  structure  is  different.  The  Shvab-Zeldovich 
relation  (13)  shows  that 


(38) 


Z.  = 


(Tg-T)  -  (2-L)Z 

to  00  G 


T  -T 
a>  a 


-1/2 

no  longer  vanishes.  As  a  consequence  the  flame  sheet  now  has  thickness  D# 
and  the  structure  is  governed  by  a  linear  equation  that  can  easily  be  integrated. 
We  find 


(39) 


,  L(l+Z  ) 

M  =  Ma[l  +  (D»Yq,)  XF^T-y 


+  o(d;1/2)] 


and  the  structure 


T  -  +  1  xir5-  v  ■><!-*  ^‘5> 1  -  »(i>;l/2). 


(40)  *  ~s  ’  "*  “0  ‘  l+z 

l/o 

where  now  5  =  D#  (r-l). 

Restrictions  on  the  parameters  comes  from  Mg  and  Z#.  The  (small) 
mass  fraction  of  fuel  is  positive  only  for  M  >  0,  so  that  we  must  have 


(41) 


T  <  T  according  as  L  *  2; 

*  8  ? 


the  positivity  of  Z*  imposes  the  additional  restriction 
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(42)  V  Ts  <  (L-2^ 

Two  angular  regions  are  thereby  determined  in  the  parameter  plane  (Fig.  3);  one 
coincides  with  the  region  where  the  Burke-Schumann  solution  exists  but  the  weak- 
burning  solution  does  not,  and  the  other  with  the  region  where  the  reverse  is  true. 

3.  The  Nearly  Adiabatic  Flame  for  Le^  ^  1. 

Having  determined  the  possible  ends  of  the  M,D-response  curve,  we  now  turn 
to  its  shape  in  between.  Here  activation-energy  asymptotics  are  particularly  useful., 
though  a  general  discussion  is  still  quite  complicated  (see  Law  1975)  requiring^ 
extensive  numerical  calculations.  For  the  moment  we  shall  restrict  attention  to 
the  special  case  considered  by  Buckmaster  (1975)  in  which 

(43)  Ta  “  T»  =  k/e  as  0  -*  «  , 

where  k  =  0(l).  A  great  deal  of  qualitative  insight  into  the  nature  of  diffusion 

flames  can  thereby  be  obtained  with  a  minimum  of  computations.  The  most  important 

aspects  of  the  general  case  will  be  dealt  with  later,  albeit  for  e-*»  still. 

It'  is  convenient  to  fix  everything  but  T  such  that  its  limiting  value 

s 

T  +  L-2  satisfies 
00 

(44)  T  <  min{T  ,T  +  (2-L)Z  } 

S  00  00  00 

to  ensure  the  existence  of  both  the  frozen  and  Burke-Schumann  limits.  Restriction 

to.  small  values  of  T  -T  means  that  in  both  limits  the  combustion  field  is  nearly 

a  oo 

adiabatic,  and  in  fact  that  will  be  true  along  the  whole  of  the  response  curve. 

The  unsteady  formulation  is  reinstated,  with  Lewis  numbers  again  arbitrary. 

Global  Shvab-Zeldovich  relations  do  not  exist  and  it  is  necessary  to  re-examine 
the  complete  system  (1-3). 


VI. 


For  slow  variations,  the  only  information  needed  from  such  a  system  in  Sec.  II 
was  the  change  in  enthalpy  T  +  Y  up  to  the  flame  sheet.  Here  the  changes  in 
hoth  T+2Y  and  T+2Z  are  needed,  tout  the  procedure  is  the  same  except  for 
integrating  from  the  supply  (as  in  Sec.  IV*)  and  taking  account  of  the  areal 


factor  r 2 . 

We  find 

(45) 

K(Ta-T„-2Y„)*y(Y)  -  -g)», 

(46) 

M(  T.-M.-S.)  =^<Z)  -  rj( f  +  §. 

to  0(1/0),  where 

0$(Y)  **r2  -r^[p(2Y-T#-2Y*)  ]dr 

1  9T 

is  to  toe  evaluated  to  q(1)  only;  stars  refer  to  values  just  behind  the  flame 
sheet.  ^(Z)  will  not  toe  evaluated  since  the  corresponding  equation  is  only 
needed  to  0  (l ) . 

So  far  as  leading  terms  are  concerned,  we  may  concentrate  on  the  steady 
version  of  equations  (2),  (3a)and  (3b)  where  pv  is  the  function  (10)with  M 
now  dependent  on  x  .  Beyond  the  flame  sheet  the  temperature  must  be  constant 
at  T  .  Suppose,  on  the  contrary , that  it  increases  from  the  flame  to  infinity. 
Ifcen  heat  from  the  infinite  resevoir  would  maintain  the  flame  temperature  above 

T  =  T  ,  which  is  a  contradiction.  The  combustion  field  is  therefore  similar  to 

a  ® 

that  for  the  plane  deflagration  wave  treated  in  Ch.  II.  The  temperature  increases 

from  T  at  the  supply  to  T  at  the  flame,  beyond  which  it  stays  constant. 

0  00 

The  reaction  in  1  <  r  <  r*  is  frozen,  while  there  is  equilibrium  with 
(47)  Y  =  0  in  r„  <  r  <  «. 

In  fact  the  latter  holds  to  all  orders. 
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With  this  picture  in  mind  the  solution  in  the  frozen  region  is  easily  derived. 
We  find 

^M(l/r#-l/r)  *»rtM(l/r#-l/r) 


(48) 


where 


(49) 


T  =  Ts-L+LeM^1_1/r^+o(l).Y=l-e' 


+o(l)  Z=Z#e 


+o(l)  , 


r*  =  M/(M+to  L/2) 


These  satisfy  the  supply  conditions  and  give  T  =  T  (to  order  1) ,  Y  -  0  at 

QO 

r  =  r#.  For  consistency  r*  must  he  greater  than  1,  i.e. 

< 

L  <  2 

is  required:  the  heat  conducted  hack  into  the  supply  sphere  must  he  less  than 
that  released  in  the  flame.  Similarly 

-M*q/r* 


(50) 


(51) 


=  ♦  Z.(e-H’n'r-e'M,"/r*)] /(1-e 


)+o(l). 


T  •  T  +  C(l-e"M/r)/0  +  o(l/e) 


in  the  equilibrium  region,  because  chemistry- free 


equations  are  still  satisfied.  (The  constant  T  does  not  count  as  a  leading 
term.)  The  unknown  constants  Z#  and  C  can  he  found  (in  terms  of  M)  from 
the  basic  equations  (45)  and  (46)  of  slow  variations,  which  yield 


(52) 


where 


(53) 


C  =  K  -bM^ dM/ dx , Z*  =  (l+Z  )e  ^  *-l 


b  =  i 


£°(2/L)  2(Le£,/2)^  - 


^-L  +  Le 


Lei  _d 
d£ 


d£ 


W 


The  chemical  reaction  determines  the  "burning  rate  at  a  given  temperature 
so  that  we  may  look  to  the  structure  of  the  flame  sheet  for  another  relation 
between  M  and  the  perturbation  of  the  flame  temperature,  represented  by  C. 
Writing  as  usual 

r  ■  r#  +  5/6  ,  T  -  T..+  t/0  +  o(l/6), 

(54) 

Y  =  y/e  +  o(1/q) ,Z=  Z#  +  z/6  +  o(i/e). 


leads  to  the  equations 


(55) 


where 


(56) 


-32t/3C2 


t/T£ 


Dy  e 


-2  -on. 

D  *  D2»e  e  , 


while  matching  on  either  side  of  the  flame  gives  the  boundary  conditions 

(57)  t  »  2M5/r2  +  0(1),  y  =  -fe/r2  +  0(1)  as  , 

-M/r# 

t  =  C(l-e  )  +  o(l),  y  =  o(l)  as  5^  +  ». 


It  follows  that  there  is  a  local  Shvab-Zeldovich  relation 

(59) 


-M/r# 

t  +  2y/£  =  C  (1-e  ), 


which  enables  the  equation  for  t  to  be  integrated  to  give 

(60) 


-M/r#  0  9/T  k  . 

C(l-e  )  =  T2jUi(2ir0  e  °7D£z*r*T„  ). 


Elimination  of  C  between  this  and  the  result  (52a)  yields  an  equation  for 


M,  namely 


VI.  15 


_  -M/r#  -1  ?  p  o  9/T  .  j. 

(61)  bw  dM/dt  =  k-(l-e  )  r to(2M  6  e  /D;£.  Z#r*T* ) 

'Where  Z#  and  r#  are  given  by  equations  (52b)  and  (49). 

We  shall  first  describe  the  steady  state,  determined  by  setting  the  right- 
hand  side  of  this  last  equation  equal  to  zero,  and  then  discuss  its  stability. 
Different  values  of  k  lead  to  different  steady-state  responses,  i.e.  curves 
of  M  versus  D  with  all  other  parameters  fixed  (Fig.  4),  though  all  join  the 
frozen  limit  (17)  to  the  Burke-Schumann  limit  (25) .  In  between,  a  curve  is 
monotonic  or  S-shaped  depending  on  the  value  of  k.  When  there  is  a  gain  of  heat 
from  infinity  (k  <0)  or  when  the  loss  is  sufficiently  small,  the  response  is  t 
monotonic.  But  for  sufficiently  large  loss  the  burning  rate  is  triple-valued  over 
some  range  of  Damkohler  number.  Conditions  under  which  heat  is  provided  to  the 
environment  sure  particularly  important,  so  that  the  S-shaped  response  is  of 
special  interest. 

Explicit  formulas  can  be  given  for  the  turning  points  on  the  S  when  k  is 
large.  At  the  lower  point 

0  ji p 

(62)  M~*n(2/L)  +  UT2/k,  D  ~2(4/e) 4l£e2e  7(£n2/L)2  ^Z*k\ 

while  at  the  upper  point 

2  Me 

(63)  M  -  Mfi  -  T~Le  ®/2k, 

0  jC£ 

D  ~4(M  +£.n  L/2)Ve  “k  exp  {-k[  l-(l+Z  )"**  ]/T2}  /LfrnM?T6e  *, 

c  00  “  oo 

where 

(64)  M  =  —  £n(l  +  Zj  +  An(2/L) 

e  -ni 

is  the  Burke-Schumann  value  (25)  corrected  forwj#  1. 


l6 


Conclusions  about  the  instability  of  the  steady  state  can  easily  be  drawn 
from  the  full  equation  (61) ,  which  holds  in  the  strip  M  <  M  <  M  where  Z#  >  0 
and  1  <  r*  <  ».  The  right-hand  side  is  positive  to  the  right  of  the  steady-state 
response  curve  and  negative  to  the  left  (Fig.  5),  while  b  £  0  according  as 

1.  Consider  first  £  <1.  If  the  burning  state  lies  off  the  steady-state 

S' 

curve,  then  M  will  change  in  the  direction  of  the  arrows  in  Fig.^so  that  the  middle 
branch  of  the  S  is  unstable.  Now  suppose£>l;  the  arrows  must  be  reversed  so 
that  then  the  upper  and  lower  branches  of  the  S,  and  the  whole  of  the 
monotonic  curve,  are  seen  to  be  unstable. 

Nothing  may  be  concluded  about  stability  however,  since  only  a  special  class 
of  disturbances  is  being  considered.  Indeed,  it  is  generally  believed  that  the 
middle  branch  of  the  S,  which  is  stable  to  the  present  disturbances,  is  in 
fact  unstable.  Whether  that  is  true  or  not,  certainly  the  top  and  bottom  branches 
are  more  important  than  the  middle  branch  because  for  sufficiently  large  and 
small  values  of  D  they  represent  the  only  possible  steady  states.  The 
instability  predicted  for  >1  is  therefore  the  most  interesting  result,  and  it 
is  analogous  to  the  instability  of  the  premixed  plane  flame  (Ch.  II )  to  one¬ 
dimensional  disturbances  when  the  Lewis  number  is  greater  than  one.  Plane  flames 
are  also  unstable  when  the  Lewis  number  is  less  than  one  if  three-dimensional 
disturbances  whose  wavelengths  are  large  compared  to  the  flame  thickness  are 
permitted  (Ch.  V),  but  this  result  clearly  has  no  analog  for  a  spherical  flame. 

If  it  is  assumed  that  the  stability  results  deduced  here  for  o£<l  are 
generally  true,  i.e.  when  all  possible  disturbances  are  accounted  for,  the 
middle  branch  is  unstable  and  the  other  two  are  stable,  then  the  turning  points 
of  the  S  are  ignition  and  extinction  points.  The  asymptotic  results  (62)  and 
(63)  then  give  ignition  and  extinction  conditions,  respectively.  In  any  event, 
there  Is  no  reason  to  doubt  their  validity  for  1. 


General  Ignition  and  Extinction  Analyses  with 


d  = 


We  proceed  now  to  the  general  steady  state,  i.e.  arbitrary  values  of  T  -  T  , 

8.  oo 

with  emphasis  on  the  top  and  bottom  branches  of  S-responses.  As  in  Sec.  2  the 
qualitative  nature  of  the  solution  is  not  affected  by  the  Lewis  numbers  so  that, 
since  the  analysis  is  complicated  enough  forX  -  =  1,  we  shall  restrict 

attention  to  these  values;  there  is  anyway  a  lack  of  stability  forjC  >  1  and  some 
question  of  it  forj£<  1,  at  least  for  T»  close  to  T&  (Sec.  3). 

We  have  to  deal  with  the  Shvab-Zeldovich  relations  (12)4(i3)  and  the  temperature 
equation  (15).  The  requirement  (44)  is  still  imposed,  for  the  same  reason.  The 
additional  requirement 


(65)  T^  Tft 

ensures  an  S-response,  as  is  suggested  by  Sec.  3  where  such  a  response  was  only 
obtained  for  sufficiently  large  k  =  8(T  -T  )>  0.  fwhen  T  is  less  than  T  , 

£L  oo  "  od  8. 

the  second  part  of  the  original  requirement  (44)  is  actually  implied  by  the  first J 
The  lower  branch  will  be  considered  first,  the  main  goal  being  to  locate 
the  ignition  point,  which  is  a  fundamental  characteristic  of  the  burning  response. 
Sec.  3  suggests  that  the  burning  will  take  place  far  from  the  supply  on  the  whole 
lower  branch,  i.e.  with  M  close  to  its  frozen  value:  the  ignition  point  (62) 
has  this  property,  albeit  the  measure  of  closeness  is  l/k  rather  than  the  1/9 
we  must  expect  now. 

At  any  finite  value  of  r  the  combustion  field  is  frozen  to  all  orders,  so 

we  may  write 


(66) 


T  *  T  -L  + 


M( 1-1/r) 
Le  + 


M  (1-1/r) 

9  -LM1L(l-l/r)e  h 


o(9_1) 


where 

M  *  Mw  +9"1M1  +  o( 9_1) • 


(67) 
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T  and  T'  take  the  correct  values  at  r  =  1  but  the  condition  T  -»-T  as 

00 

r  -*■  00  is  violated,  shoving  that  a  new  coordinate  is  required  there.  It  turns  out 
to  be 

(68)  R  =  r/0 

and  then  the  expansion 


(69)  T  =  +  6_1t(R)  +  o(0-1) 

leads  to  the  structure  equation 

,  .  „  ..  (T  -T  )M  t/T2  -  -0/T  /2 

(70)  “.o.-oz.e-H  » 

n 


subject  to  the  boundary  conditions 


M 


(71) 


t  =  Le  W(M--M, ^R)  +  o(l)  as  R  0,  t  =  o(l)  as  R  ->  », 


the  first  of  which  comes  from  matching  with  the  frozen  expansion.  With  all 
parameters  specified,  including  Dw  the  perturbation  t  can  be  found  (numerically) 
without  using  the  M^-term.  In  that  way  can  be  determined  as  a  function  of 

the  DamkShler  number. 

By  suitably  scaling  D  and  M-.  the  relation  between  them  can  be  made  to 

W  ± 

depend  on  a  single  par amt er 

M 

(72)  6  =  (VTJ/Le  W* 


The  scaling  is 

2  M  ? 

(73)  =  a  Dv,  Nw  =  aM1/Mw  where  a  =  \Le  W/T,, 


and  then  the  new  variables 
(74)  u  =  Pt/T^  ,  p  =  R/a 
lead  to  the  canonical  problem 
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(75)  d^u/dp^  =  ^(u-g)eU^p, 

(76)  u  =  N^-l  +  o(l)  as  p->- 0,  u  =  o(l)  as  p  «. 


It  is  found,  albeit  numerically,  that  the  graphs  of  Nw  versus  A  v  display  the 
ignition  phenomenon  if  and  only. if  fj  is  positive.  Curves  for  various  values  of 
6  are  shown  in  Fig.  6;  note  that  they  are  independent  of  Zm  . 

Hie  burning  rate  at  the  extinction  point  (63)  of  the  nearly  adiabatic 
flame  lies  close  to  the  Burke-Schumann  value  when  k  is  large.  This  suggests 
that  the  general  extinction  curve  can  also  be  uncovered  by  an  analysis  which 
considers  o(0  )  perturbations  of  the  burning  rate,  now  from  its  Burke- 

Schumann  value. 

The  flame  sheet  has  the  location  (26)  while  on  both  sides  of  it  the  reaction 
is  frozen  to  all  orders.  Thus 


(77)  T  = 


T~  =  Tg-L+Le 


“ed-i/r)  .1, 


+  0  M^Ld-l/r)  e 


M  (l-l/r) 


+  o( 0_1)  for  r  <  T 


+  _Me/r 

T  =  Ta+(Too-Ta)e  +o(l)  for  r  >  r*. 


where 


(78) 


M  =  Me  +e-1M1  +  o(  0’1). 


Only  is  unknown  and  for  that  we  must  turn  to  the  structure  of  the  flame 

sheet.  With 


(79)  r  =  r*  +  5/0 


the  expansion 


(80)  T  =T,  +e-1t(  5)  +  o(  0_1), 


where  T#  is  the  flame  temperature  (27),  leads  to  the  structure  equation 
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it  T  -T  M  5  M.  T  -T  Mr  M.  t/TT  .  -e/T,/4 

(8»  ^  .  -Det  H  t  ♦  (3V-  2)(-|*  .  jl,  ]  .  *.  *-  D9-3e 

subject  to  the  boundary  conditions 

T  -T  M  £  > 

t  «  (2  -  3^-^  [-J-  +  (1  +  °d)  as  €-*■-«» 

<=»  r-  * 

(82) 

..  T  -T  M 

dlj  =  ~(i?Z?~2  +  o(l)  as  €■*■  +0°  » 

r* 

vhich  comes  from  matching.  This  should  be  compared  with  the  problem  (33),  (34) 

and  (  3fi .  the  only  difference  being  that  here  the  Arrhenius  factor  is  not  approximated 

t 

by  a  constant  since  the  perturbations  are  as  large  as  0(e_1). 

Again  the  problem  determines  and,  by  suitable  scaling,  the  relation  between 

it  and  Dfi  can  be  made  to  depend  on  the  single  parameter 

(83)  Y  -  1  -  (Ta-To)/d+z«o)  =l-g. 


Ae  =  ri^De/^,  Ne  =  (1+y)Hj/tJ, 


u  =  t/T*  -yp  ,  p  =  (Hgg/rJ^/r,)/^ 


leads  to  the  canonical  problem 


(86)  d2u/dp2  =  -  A  (u2-p2)exp(u  +yp). 


(87)  u  =  P+  Ng  +  o(l)  as  p-H-oa  ,  du/dp  =  -1  +  o(l)  as  J>L*  +  oo 

It  is  found  (though  again  numerically)  that  the  graphs  of  Ng  versus  display 

the  extinction  phenomenon  if  and  only  if  y  is  less  than  1,  i.e.  g=  1  -v  is 
positive.  Curves  for  various  values  of  Y  are  shown  in  Fig.  7.  Note  that  the 

3  e/T# 

actual  Damkohler  number  for  extinction  is  0(eJe  )  as  compared  to 

2  0/T« 

0(8  e  )  for  ignition,  so  that  (since  T*  >  T, )  the  extinction  point  lies  to 
the  left  of  the  ignition  point,  as  it  should. 
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It  is  widely  believed  that  the  middle  branch,  where  M  is  a  decreasing 
function  of  D,  is  not  of  practical  interest  because  the  burning  states  there  are 
unstable.  However,  there  is  no  hard  theoretical  evidence  for  this  other  than  the 
slow-time  instability  for  <  1  described  in  Sec.  3.  Certainly  the  notion  that 
for  a  laboratory  flame  an  increase  in  pressure  could  cause  a  decrease  in  burning 
rate  is  intuitively  difficult  to  accept  but,  be  that  as  it  may,  there  are  at 
least  two  good  reasons  for  studying  the  middle  branch.  Any  definitive  answer 
to  the  stability  question  can  surely  come  only  after  the  nature  of  the  steady  i^fcate 
has  been  determined;  and  anyway  that  nature  is  sufficiently  different  from  what 
we  have  encountered  before  that  some  account  of  it  should  be  given.  Nevertheless, 
only  a  broad  sketch  will  be  presented. 

.  It  is  clear  from  the  ignition  and  extinction  analyses  that  the  flame 
temperature  on  the  middle  branch  must  span  the  range  T  to  T  +  Z  (T  -T  )/(l+Z 

CO  00  00  3.  00  oc 

and,  since  it  is  the  maximum  temperature  at  both  ends,  we  may  expect  it  to  be  so 
ix  between.  On  either  side  of  the  flame  sheet  the  reaction  is  frozen,  just  as  in 
the  extinction  analysis.  The  key  difference  is  that  here  Y#  and  Z#  do  not 
both  vanish;  we  shall  first  suppose  that  neither  does. 

It  is  more  convenient  now  to  specify  M  and  determine 


C88)  D  -  DQ[l  +  O(0_1) ] 


as  part  of  the  solution.  The  expansions  (77 )  are  therefore  replaced  by 


(39)  T 


T"  *  T  -L+LeM(l_1/r^  +  o( 1 )  for  r  <  r#. 


T  =  T  +  Le^(l-e”M//r)  +  o(l)  for  r  >  r, , 

Ob 
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if  we  anticipate  a  property  of  the  flame  structure,  namely  that  the  heat  generated 
by  the  reaction  is  conducted  equally  to  the  two  sides.  Continuity  of  temperature 
to  leading  order  at  the  flame  sheet  now  determines 

M  -M 

(90)  r*  =  M/to  2/(l+e  W  ) 
and 

(91)  T*  =  (Ts+Too-L+LeM)/2, 


which  covers  its  range  as  M  increases  from  Mw  to  M  . 

The  sturcture  is  investigated,  in  the  usual  way,  by  means  of  the  trans format irr. 
(79)  and  the  expansion  (80  )  and  we  find  ^ 


(92) 


$  -  -Vt/T|. 


m 


"  o0Y.z*9‘le'e/T* 


where  Y#  and  Z#  are  determined  by  the  Shvab-Zeldovich  relations  (12)  and  (13) . 
Integration  yields 

P  o  t/T* 

(93)  (f|) 2  =  C  -  2DmT2e  , 


where  C  is  a  constant;  from  which  it  follows  that  the  temperature  gradients  on 
the  two  sides  of  the  flame  sheet  are  of  equal  magnitude  but  of  opposite  sign,  as 
anticipated  earlier.  [Matching  shows  that  C  =  (Le^M/r#)exp(-M/r#) ,  but  we 
shall  not  need  this  result.] 

Another  integration  gives  only  a  relation  between  the  maximum  value  of  t 
and  its  location;  to  determine  it  is  necessary  to  take  the  expansion  (80) 

to  one  more  term.  No  useful  purpose  is  served  by  presenting  these  rather 
complicated  algebraic  details.  They  are  described  by  Lindn  (197*0  for  the 
counterflow  flame  and  Law  (1975)  has  adapted  Linan's  analysis  to  the  present 
problem.  Another  reason  for  not  pursuing  the  determination  of  is  that 


is  only  affected  to  0(l)  by  it,  T#  being  a  known  function  (91  )  of  M. 

Indeed  the  expected  property 

dDg/dM  <0  as  8  +  « 

follows. irrespective  of  the  value  of  D  .  since 
7  m 

(95  )  dT*/dM  >  0. 

The  solution  derived  here  is  valid  only  so  long  as  0  <  Y^Z*  <  1  and 
r#  >1.  These  inequalities  are  satisfied  automatically  with  the  exception  of  * 

Yp,#>  0,  which  requires 

(96)  M  <  £n(T  -T  )/L. 

a  °° 

Since  M  must  be  greater  them  Mw  the  condition  for  there  to  be  a  range  of  validity 
is 

(97)  T«  <  Ta  "  ^ 

the  range  will  extend  up  to  Mg  if 

(98)  <  V(1+ZJ* 

These  possibilities  are  shown  in  Fig.  3. 

When  this  last  condition  is  not  satisfied,  a  question  arises  of  completing 

the  range  Mw  <  M  <  Mg.  The  clue  is  the  vanishing  of  Y*  at  the  end  of  the  range: 

for  the  remaining  values  of  M  (and  all  of  them  when  T  -T  <  1)  the  mass  fraction 

a  <*>  ' 

of  fuel  is  oCe'1)  at  the  flame  sheet.  This  leads  to  a  similar  analysis,  albeit 
with  a  different  flame  structure  (cf.  Sec.  6),  which  completes  the  middle  branch. 

On  it  Z#  decreases  from  at  one  end  to  zero  at  the  other  so  as  to  join  the 


VI. 


ignition  and  extinction  branches.  In  that  connection  it  should  be  noted  that 
the  middle  branch  does  not  join  directly  to  the  extinction  branch  when  the  original 
solution  completes  it.  Although  Z#  vanishes  at  M  ,  Y#  does  not;  so  that 
a  transition  flame,  for  which  Y#  decreases  to  zero,  is  needed.  All  these 
details  are  described  by  Linan  (197*0  for  the  counterflow  problem. 

6.  The  Monotonic  and  Other  Responses. 

The  last  two  sections have  been  concerned  with  the  S-shaped  response  which 
occurs  when  the  combustion  field  supplies  heat  to  the  ambient  atmosphere.  When 

(99  )  T  >  T 

oo  a 

the  flux  is  in  the  opposite  direction  and  we  may  expect  the  response  to  be 
monotonic  (if  only  on  the  basis  of  Sec.  3).  Fuel  drops,  to  which  the  theory  will 
be  adapted  in  Sec.  8,  are  usually  burnt  to  provide  heat  to  their  surroundings 
so  that  from  the  point  of  view  of  that  precise  application  the  present  discussion 
is  not  important. 

The  temperature  increases  again  beyond  the  flame  sheet  to  give  an  equilibrium 
region,  where 

(100)  Y  =  0  for  r  >  r# 

must  hold  to  all  orders.  It  follows  that  the  expression  (89a)  is  still  correct, 
but  that  the  expansion  (89b)  must  be  replaced  by  the  asymptotically  exact  result 

(101)  T  =  T+  =  T  +  (T  -T  )(l-e_M/r)  for  r  >  r*. 

Continuity  at  the  flame  sheet  then  determines 

r,  =  M/ In  t(Ta-T+  LeM)/2  ] 

T„  =  Ta  +  2(Too-Ta)/(LeM-Too+Ta). 


(102) 


The  monotonicity  of  the  response  is  now  clear  without  going  into  details  of  the 
structure,  by  extracting  just 

(103)  to  D0  =  e/T*  +  2  in0  t  0(1) 
in  analogy  with  the  expansion  (36).  Since 

(104)  dT*/dM  <  0 

under  the  condition  (99  ) ,  the  leading  term  Dq  in  the  DamkCfhler  number  is  an 
increasing  function  as  required. 

As  M  increases  from  its  frozen  value,  r*  decreases  from  <»  and  there  ^ire 
two  possibilities.  If  the  inequality  (36)  is  satisfied  (see  Fig.  3),  r*  reaches 
the  value  (26  )  as  M  tends  to  Mg;  the  combustion  field  is  then  identical  to 
the  Burke-Schumann,  in  particular  the  oxidant  no  longer  penetrates  through  the 
flame.  Our  solution  joins  the  frozen  limit  to  the  Burke-Schumann  limit,  correct 
descriptions  of  the  approaches  to  which  follow  the  same  lines  as  Sec.  4. 

The  second  possibility  arises  when 

(105)  0  <  T  -  T  <  ( L-2 )  Z 

00  S  op 

(see  Fig.  3).  Then  M  is  smaller  than  M  and  as  M  reaches  it  the  flame 
settles  down  on  the  surface  of  the  supply  sphere.  Our  solution  now  joins  the 
frozen  limit  to  Buckmaster's  limit,  though  we  shall  not  describe  the  approaches  here. 

The  discussion  so  far  in  this  section  can  be  justified  by  the  fuel  drop,  albeit 
under  conditions  of  little  practical  interest.  The  same  is  true  for 

(106)  To#<  Ts, 

though  the  inequality  can  no  longer  hold  over  the  whole  response  of  the  drop  (see 
the  next  section) .  For  this  reason,  no  analysis  for  the  lower  half  of  Fig.  3  has 
been  published  although  it  is  not  difficult  to  see  how  it  would  go. 
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While  there  is  no  frozen  limit  now,  it  turns  out  that  both  equilibrium 
limits  exist,  provided 


(107)  T„-  Ts>  (L-2)Z„. 

Otherwise,  no  limits  exist  and  we  must  conclude  that  there  is  no  solution.  Indeed, 
application  of  the  Shvab-Zeldovich  relations  (12)  and  (13)  at  r  =  1  shows  that 
the  fuel  and  oxidant  fractions  there  would  be  negative  for  L  >  2  •  and  for 
T^-T^  <  -2(1+2*,),  respectively,  even  for  finite  0  ,  which  leaves  only  a  small 

triangle  of  doubt. 

t 

It  is  reasonable  to  suppose  that,  when  the  conditions  (106)  and  (107)  are  satisfied, 

Ihe  response  takes  the  form  of  a  C  joining  Buckmaster's  limit  to  the  Burke-Schumann 

limit  (which  is  always  higher).  Indeed,  the  analysis  sketched  in  Sec.  5  will 

still  apply;  only  the  conditions  under  which  the  basic  solution  must  be  modified 

will  change:  r#  >  1  is  no  longer  automatic  and  the  basic  solution  applies  over 

at  least  part  of  the  range  M_,  <  M  <  Jt  only  when  L  <  1. 

s  e 

The  leftmost  point  of  the  C-shaped  curve  represents  both  ignition  and 
extinction  conditions.  Since  the  response  determined  by  Sec.  5  has  negative  slope, 
the  burning  rate  there  will  again  lie  within  0(l/9)  of  the  Burke-Schumann 
value  Me-  Hence  the  extinction  analysis  of  Sec.  h  locates  these  conditions. 

7.  The  Burning  Fuel  Drop. 


In  our  discussion  of  the  spherical  diffusion  flame  we  have  supposed  the  seven 
parameters 

*,.L,Js.Ks,T0.,Y.,^ 
to  be  given;  indeed,  we  have  taken 


Js  »  1,  Ks  =  0,  Ym  =  0 
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for  simplicity.  If  L  is  the  latent  heat  of  the  vaporization  of  the  fuel  (taken 

to  be  constant  for  the  temperature  range  of  interest),  then  these  are  in  fact 

prescribed  in  fuel -drop  burning  with  the  exception  of  T  ,  in  place  of  which  the 

s 

Clausius-Clapeyron  relation 

(108)  Ys  *  (Ts/T^.6  exp  eU/^-l/Tg) 

is  substituted.  Here  is  the  boiling  temperature  corresponding  to  pc  i.e., 

in  the  notation  of  equation  (III.16), 

(109)  exp(-8/Tb)  =  pc/k.  f 

Now  Y  may  be  calculated  in  terms  of  the  existing  parameters  and  M  by  means 
of  the  Shvab-Zeldovich  relation  (12)  ,  so  that  it  may  be  eliminated  to  give 

(110)  2(T  /T.  )eexp6(l/T.-l/T  )  =  2-L+(T  -T  +L-2)e~M 

SO  D  S  00  S 

as  the  equation  determining  T  .  The  new  feature  therefore  is  that  T  is  no 

s  s 

longer  fixed  but  must  be  calculated  at  each  point  of  the  response  curve.  Once 
Ta  is  found,  the  parameter  plane  of  Fig.  3  may  be  used  to  find  the  appropriate 
structure  for  determining  D. 

The  properties  of  the  resulting  responses  are  not  obvious  and  it  helps 
to  characterize  the  loci  in  the  parameter  plane .  First  note  that  there  is  always 
a  weak  burning  limit,  whatever  the  values  of  L,T  ,Z  ,0  and  T.  .  That  is,  the 

oo  00  D 

equation  (110)  invariably  has  a  solution  T  =  T  (say)  when  M  is  set  equal  to 

s  w 

in  fact,  T^  <  Tm  as  we  should  expect  in  the  absence  of  weak  burning 
for  T  >  T  .  Next  note  that  there  is  always  a  Burke-Schumann  limit 
T  *  T  (say)when 

O  ™ 


(111) 


T  >  (L-2)Z  ; 

w  00 
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that  T  <  T  +(2-L)Z  ,  corresponding  to  a  point  above  the  quilibria  divider; 

0  00  00 

and  that  Te  J  according  as  the  point  is  above  or  below  the  adiabatic  line. 

These  results  follow  on  setting  M  =  Mg  in  equation  (110).  Finally  note  that. 

when  M  =  M  ,  the  equation  has  the  solution  T  =  0.  While  this  is  not  relevant 
s  s 

when  the  inequality  (1+4)  holds,  it  provides  a  termination  point  when 

(112)  <(L-2)Zm. 

The  possible  loci  are  illustrated  in  Fig.  8,  where  dashes  indicate  that  the 
corresponding  lines  may  or  may  not  be  crossed  (depending  on  T^  and  0).  The  t 
arrows  show  the  direction  of  increasing  M  and  follow  from  the  property 

(113)  dT  /dM  S  0  according  as  T  I  T  +  L-2 

s  s 

of  equation  (110).  I  The  condition  T  <_  T.  ,  i.e.  Y  <  1,  is  always  satisfied.] 

Our  knowledge  of  the  response  curves  for  fixed  points  in  the  paramter  plane  now 
suggests  that  the  location  of  the  point  w  determines  whether  the  fuel-drop 
response  is  S-shaped  or  monotonic.  Proof  lies  in  the  inequalities  (95) &  (104), 
which  continue  to  hold  when  T  varies  with  M  according  to  the  equation  (110) . 

In  short,  qualitative  results  for  fuel-drop  burning  can  be  deduced  from  the 
analysis  of  general  spherical  flames  given  in  previous  sections.  However, 
quantitative  results  such  as  the  graphs  in  Figs .  6  &  7  cannot  be  carried  over : 
the  variations  of  T  ,  while  they  are  only  o(l/9)  on  the  upper  and  lower  branches, 
will  modify  the  ignition  and  extinction  values  (Ludford  &  Hormandia  1978). 

These  remarks  apply  to  steady  burning;  the  unsteady  analysis  of  Sec.  3  is 
not  directly  applicable  since  no  account  is  taken  of  conditions  in  the  interior 
of  the  drop.  Some  of  the  heat  flux  at  the  surface  now  accounts  for  changes  in  the 
interior  temperature,  so  that  there  is  a  connection  between  it  and  the  surface 
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temperature .  As  a  consequence,  L  is  the  sum  of  the  latent  heat  of  vaporization 
and  a  complicated  integral  term  depending  on  T  .  Such  difficulties  do  not  arise 
for  the  counterflow  flame,  but  there  the  serious  compromises  with  the  fluid 
mechanics  which  are  necessary  cast  doubt  on  the  results. 
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